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Multiscale modeling at small length scales (10 9 to 10" 3 m) is discussed for aluminum 
matrices with embedded particles. A configuration containing one particle surrounded by 
about 50 grains and subjected to uniform tension and lateral constraint is considered. The 
analyses are performed to better understand the effects of material configuration on the 
initiation and progression of debonding of the particles from the surrounding aluminum 
matrix. Configurational parameters considered include particle aspect ratio and orientation 
within the surrounding matrix. Both configurational parameters are shown to have a 
significant effect on the behavior of the materials as a whole. For elliptical particles with the 
major axis perpendicular to the direction of loading, a particle with a 1:1 aspect ratio 
completely debonds from the surrounding matrix at higher loads than particles with higher 
aspect ratios. As the particle major axis is aligned with the direction of the applied load, 
increasing amounts of load are required to completely debond the particles. 


I. Introduction 

M ULTISCALE modeling provides an efficient means of interrogating deformation and fracture of metallic 
materials from the micro- to the nano-scales. Cohesive zone models are an important technique used in 
multiscale modeling of grain-scale events and can be used to study grain boundary decohesion of metallic materials 
under loading. Similar grain-scale simulations can be used to quantify the 
effect of particles on load transfer and fracture initiation within the 
metallic materials. 

Determination of the interaction (particularly, the load transfer) 
between particles and a surrounding metallic microstructure (Figure 1) is 
a topic of considerable attention within the materials processing 
community. Deformation and fracture at length scales such as those 
associated with the particle may not be adequately modeled by traditional 
continuum mechanics-based concepts. Rather, the atomic-level events 
related to bond breakage, dislocation motion, and void nucleation must be 
understood and then used to predict debonding at the grain scale. This 
process requires an understanding of material behavior at the atomic level 
as well as development of means to bridge length scales. Figure 1: Particles embedded in a 

In this paper, multiscale modeling of deformation and fracture of metallic matrix subjected to far-field 
metallic materials with embedded particles will be discussed. Cohesive loading 
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zone models are used to estimate the effect of the atomic-scale events in finite element (FE) -based micromechanics 
models of an elliptical particle embedded within an aluminum matrix. The micromechanics model is used to 
simulate the debonding along the grain boundaries of the aluminum matrix and between the grains of the aluminum 
matrix and the particles. The effects of particle size, shape, and orientation on the debonding of particles from 
aluminum grains will be determined. Additionally, the effects of the debonding of the particle on load transfer and 
fracture of the material as a whole will be discussed. 

II. Multiscale Modeling 

Length scales on the order of 10" 9 to 10" 3 m dominate the problem of particle-matrix debonding. Atomistic 
effects dominate at the lower bound of this range, while traditional continuum assumptions are valid near the higher 
bound. In the present analysis, cohesive zone models (CZM) 1 , which approximate the interfacial atomic traction- 
displacement relationship using bulk material properties, are used to model the response along the grain boundaries 
and between the particles and the grains of the aluminum matrix. Cohesive zone models assume cohesive 
interactions among the grains of a material and permit the appearance of fracture surfaces in a continuum. 2 One 
popular cohesive zone model is attributed to Tvergaard and Hutchinson 1 , where the normal and shear components of 
the traction and displacement are combined into single measures, r and X, respectively, so that the responses are 
coupled. 3 The coupled cohesive zone model (CCZM) given in reference 3 defines a traction potential, 

0 (S n ,S t ) = S c n \T(X')dX\ ( 1 ) 

X 


where X is a non-dimensional measure of the relative opening and sliding displacements ( S n and S t ) as defined by 
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(a) Hypothetical traction-displacement 
relationship from MD 



S n c and S t c are the critical values for the opening and sliding modes, respectively, and the initiation of a crack is 
assumed to occur once X reaches a value of unity. The present methodology for estimating atomic scale information 
within the grain scale models is shown in Figure 2. Figure 2(a) represents a hypothetical traction-displacement 
relationship describing the 
atomic-level events at the grain 
boundary, while Figure 2(b) is 
the estimated response 
embedded within the CCZM. 

The traction-displacement 

response in Figure 2(b) is 
embedded within cohesive zone 
elements (Figure 2(c)) to model 
the response of the grain 
boundaries. The cohesive zone 
elements are embedded between 
continuum finite elements on 
either side of the grain boundary 
as shown in Figure 2(d). 1 In two 
dimensions, the cohesive zone 
elements have zero width in the 
normal direction, i.e., the 
unprimed and primed nodes in 
Figure 2(c) are coincident. 


(b) Estimated traction-displacement 
relationship in CCZM 
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(c) Cohesive zone element 


(d) Cohesive zone elements embedded in the grain 
boundaries of the crack path within a grain structure 


Figure 2: Implementation of the cohesive zone models. 
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III. Micro-Scale Modeling 

Voronoi tessellations are used to represent the grains of the metallic material. The tessellations are formed by 
first randomly distributing nuclei and connecting these nuclei via straight lines to their nearest neighbors, forming a 
Delaunay triangulation. 3 4 These lines are then perpendicularly bisected to create the edges of a polygon. The 
resulting polygons define the Voronoi regions, the areas within which all points are closer to the nucleus of that 
polygon than to any other polygon nucleus. 5 In the current work, Voronoi polygons represent the individual grains 
of a metallic material surrounding the domain of the particle (see Figure 3). Particles in this study are 1/3 to 1 times 
the size of the typical matrix grain size. The algorithm by which these models are generated is outlined below: 

1. An elliptical particle and buffer zone are defined. Voronoi 
nuclei are randomly chosen such that they do not lie within 
the buffer zone, and a tessellation is formed (see Figure 4). 

2. The Voronoi boundaries that intersect the buffer zone are 
clipped to the buffer, and the clipped points are connected 
(see Figure 5). 

3. Each of the clipped points is connected to a major axis 1/3 
point of the elliptical particle, and then these new edges are 
clipped to the ellipse boundary. The points clipped to the 
ellipse are connected, resulting in new regions between the 
particle (which is distorted) and the existing regions. The 
curved edges of the particle are recovered by defining arc 
segments between the vertices. The results of these 
processes are shown in Figure 6. 

4. Small grains that form around the ellipse are merged with 
neighboring grains to make a more realistic model, and the 
model is meshed. A finite element mesh of the model in 
Figure 6 is shown in Figure 7. 


t t t w t t t 


$> 

l > 8 > 



/ 


Particle 


Surrounding metallic microstructure 


<# 

Particle orientation angle 

<# 

<# 


A 


A 

- K 


A 


Figure 3: Model of material microstructure 
with a prescribed elliptical-shaped particle. 


IV. Analysis 

The work in reference 6 is extended here to model load transfer and debonding of particles in an aluminum 
matrix. A Voronoi tessellation is generated around a prescribed elliptical-shaped particle as shown in Figure 3. The 
problem is assumed to be plane strain and is subjected to uniform tension and lateral constraint. The matrix is 
assumed to be isotropic elastic 7075 -T6 aluminum, and the particle (shaded in Figure 3) is assumed to be a relatively 
soft isotropic elastic material. Cohesive zone elements implementing two different CCZMs are placed along the 
grain boundaries: one model for the aluminum-aluminum interface, and one model for the particle-aluminum 
interface. The material properties used for the 7075-T6 aluminum grains, the particle, and the grain boundaries are 
given in Table 1. 

The aluminum 
matrix and particle 
properties include 
the modulus, E , and 
the Poisson’s ratio, 
v. The grain 
boundary properties 
include the peak 

stress, t p , the initial stiffness, k n , and the critical opening displacements normal and tangential to the interface, S n c 
and 8 t c , respectively. The parameters that define the CCZM are chosen such that when all of the grains, including 
the particle, are assigned homogeneous material properties and a very small displacement is applied, load is 
completely transferred across the CCZM and a homogeneous stress field is obtained. The peak traction for the 
particle-aluminum interface is assigned a lower value than the aluminum-aluminum interface to guarantee that the 
particle debonds from the aluminum matrix before crack initiation occurs in the aluminum matrix. The simulation 
sample has dimensions (/ x , l y ) along the x andy axes, where l x = l y = 1 mm. 

As the grains in the present work are all considered to be isotropic and of the same elastic properties, the grain 
boundaries should have no effect on the analyses before grain boundary decohesion initiates; results should agree 


Table 1: Domain and grain boundary properties 


Domain Properties 

Boundary Properties 


Aluminum 

Particle 


Aluminum- Aluminum 

Particle-Aluminum 

E 

72 GPa 

57 GPa 


0.5 GPa 

0.1 GPa 

V 

0.33 

0.32 

kn 

4x 10 9 GPa 

4 x 10 9 GPa 




o 
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1 pm 
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Figure 4: Voronoi diagram and overlaid 
particle and buffer zone. 



Figure 6: Voronoi diagram with curved 
edges of the particle recovered and 
underlaid distorted particle. 



Figure 5: Voronoi diagram with 
clipped points connected and overlaid 
particle. 



Figure 7: Finite element mesh of the model 
represented by Figure 6. 


Uniform, identical stress distribution 
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with those obtained by performing the 
analyses on particles embedded in an 
amorphous material. Figure 8 
presents the loading-direction stress 
distribution at 0.1% applied strain for 
a circular particle embedded in (a) 
amorphous isotropic aluminum and 
(b) polycrystalline isotropic 
aluminum. As expected, the stress 
distributions in Figures 8(a) and (b) 
are identical, and a stress 
concentration is formed by the 
particle. 

Legarth 6 studied the effects of 
particle aspect ratio and orientation on 
the onset and completion of debonding 
of rigid elliptical particles from metal 

matrices. In the analyses, the matrix material was considered to be amorphous, and both plastic isotropic and plastic 
anisotropic material models were used. Coupled cohesive zone models were used to model the particle-matrix 
interface. Particles with 0° and 90° orientations and aspect ratios of 1:1, 2:1, and 4:1 were considered. 
Conclusions from reference 6 include: 

1 . The onset of debonding is nearly unaffected by the particle aspect ratio. 

2. The complete debonding of the particle from the matrix requires higher strains as the major axis of the particle is 

aligned with the axis of load application. 


^ 99.2 


(a) Amorphous isotropic aluminum (b) Polycrystalline isotropic aluminum 

Figure 8: Loading-direction stress distribution at 0.1% far-field strain 
for a circular particle embedded in aluminum 
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In the current study, models with elliptical particles of three different aspect Table 2: Aspect ratios (AR) 
ratios (AR) and four angles of orientation (0) (given in Table 2) are considered. 

A summary of the problem parameters used in both reference 6 and the current 
work is presented in Table 3. The boundary conditions and loading for a typical 
(/*, l y ) model with a particle oriented at 0 degrees within a matrix tessellation are 
shown in Figure 3. Particles are placed in the center of the models and are 1/3 to 
1 times the size of the surrounding grains. In the current work, the particle 
minor axis is held constant, and as the aspect ratio of the particle increases, so 
does its area. 

Uniform far-field tension 
is applied in the form of 
0.01% strain increments to a 
maximum strain of 0.5%. 

Lateral constraint at x = 0 
and l x (see Figure 3) is 
assumed and will cause 
transverse stresses, cr x , in the 
models. The models are loaded until the peak traction, z p , of the particle-aluminum interface is reached and the 
elliptic particles begin to debond from the surrounding matrix. The far- field strain to initiate debonding is shown in 
Figure 9 as a function of particle orientation angle, 0, and particle aspect ratio, AR. These results indicate a small 
dependence (approx. 5%) of debond initiation on both particle aspect ratio and angle of orientation with the lateral 
constraint boundary condition. 

Loading of the elliptical particle models continues until the particles debond completely from the aluminum 
matrices. In this work, the particles are said to have completely debonded when tractions are no longer transferred 
across the particle-aluminum interface. As seen in Figure 10, increasing the angle of orientation of the 2:1 and 3:1 
particles increases the strain required to completely debond the particles. 

Table 4 presents the strain increase required to reach complete debonding from debonding onset for the three 
different aspect ratios for each of the orientation angles considered. For a 1:1 aspect ratio (i.e., a circular particle), a 
larger increase in strain is required to completely debond the particle once onset occurs than is required for the 2:1 


Table 3: Problem parameters used in reference 6 and the current work 



Reference 6 

Current Work 

Particle character 

Rigid 

Elastic 

Material properties 

Anisotropic 

Isotropic 

Particle orientation angles 

0° and 90° 

0°, 30°, 60°, and 90° 

Matrix configuration 

Amorphous 

Distinguishable grain boundaries 


and angles of orientation (0) 


AR 

0 

1:1 

0° 

2:1 

30° 

3:1 

60° 


90° 



Angle, degrees 

Figure 9: Far-field strain to initiate debonding (reak peak traction) for particles of three different aspect 
ratios embedded in an aluminum matrix. 
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Figure 10: Far-field strain to completely debond particles of three different aspect ratios from an 

aluminum matrix. 


Table 4: Strain increment required to completely 
debond particles of three different aspect ratios from 
an aluminum matrix once onset has occurred 



1:1 

2:1 

3:1 

0° 

0.200% 

0.080% 

0.060% 

30° 

0.200% 

0.090% 

0.060% 

o 

O 

0.200% 

0.140% 

0.100% 

^sO 

o 

o 

0.200% 

0.260% 

0.200% 


Strain in soft particle is larger 
than strain in aluminum 


and 3:1 aspect ratios. For non-circular elliptical 
particles, larger strain increases are required to 
completely debond the particles as the orientation angle 
aligns the particle major axis with the direction of the 
applied load. Figures 9 and 10 can then be used to 
estimate the strains at which initial and final debonding 
will occur. For example, particles in all three aspect 
ratios, when oriented at 45°, will begin to debond at 

0. 11% strain, and the 1:1, 2:1, and 3:1 particles should 
completely debond at 0.31%, 0.22%, and 0.18% strain, 
respectively. 

The details of the load transfer between the particles and the 
surrounding aluminum matrix and the effect of the particles on the 
fracture of the matrix are now considered for a configuration with an 
elliptical particle having an aspect ratio of 2:1 and orientation of 0°. As 
shown in Figure 1 1, the strain in the relatively soft particle is larger than 
the strain in the aluminum if the local grain boundary remains connected, 

1. e., the peak traction of the particle-aluminum interface ( t p = 0.1 GPa) is 
not surpassed. Once the peak traction is surpassed, the particle begins to 
debond from the surrounding aluminum (Figure 12(b)), and the load 
transfer across the boundary is reduced (Figure 12(a)). Stress 
concentrations form at the narrow edges of the elliptical particle as 
shown in Figure 12(a). As more strain is applied to the material, the 
particle continues to debond and the stresses in the material increase until 
the particle completely debonds from the aluminum. Further increasing 
the applied strain increases the stresses in the aluminum until the peak 

traction of the aluminum-aluminum interface ( t p = 0.5 GPa) is surpassed (at 0.20% strain), and the crack begins to 
grow in the aluminum grains (see Figure 13). As the crack propagates (Figure 14(b)), the stress in the material is 
redistributed as shown in Figure 14(a). 
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Figure 11: Loading-direction strain 

distribution for 0.1% far-field strain. 
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(a) Strain distribution (b) Grain boundary opening 

Figure 12: Loading-direction strain distribution and grain 

boundary opening for 0.15% far-field strain. 



(a) Stress distribution (b) Grain boundary opening 

Figure 13: Loading-direction stress distribution and grain 

boundary opening (crack initiation) at 0.20% far-field strain. 


V. Concluding Remarks 

Multiscale modeling is discussed for 
isotropic aluminum matrices with 
embedded particles. The behavior of the 
grain boundaries and the particle- 
aluminum interface is modeled using 
coupled cohesive zone models, which use 
bulk material properties to estimate an 
atomic-scale traction-displacement 

relationship. 

The effects of the particle size, shape, 
and orientation within the aluminum 
matrix on debond initiation and 
completion, load transfer, and fracture of 
the matrix are studied. For the case of 
uniform far-field tension loading and 
lateral constraint boundary conditions 
considered, the configuration of the 
particle has little effect on the load at 
which debonding initiates. However, the 
configuration of the particle has 
significant effect on the load required for 
complete debonding. For elliptical 
particles with the major axis 
perpendicular to the direction of loading, 
a particle with a 1:1 aspect ratio 
completely debonds from the surrounding 
matrix at higher loads than particles with 
higher aspect ratios. In addition, the load 
required to completely debond the particle 
increases as its major axis is aligned with 
the loading direction. Debonding of the 
particle precipitates cracking along the 
grain boundaries of the aluminum matrix. 


Stress reduction 

and redistribution cr, (MPa) 



Stress concentration 
at crack tip 


(a) Stress distribution 



Crack tip Through-crack 


(b) Grain boundary opening 


Figure 14: Loading-direction stress distribution and grain 

boundary opening (crack propagation) at 0.21% far-field strain. 
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